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Q ■ Abstract 

Using methods from effective field theory, an exact all-order expression for the Drell-Yan 
cross section at small transverse momentum is derived directly in q? space, in which all 
large logarithms are resummed. The anomalous dimensions and matching coefficients 
necessary for resummation at NNLL order are given explicitly. The precise relation be- 
tween our result and the Collins-Soper-Sterman formula is discussed, and as a by-product 

r | the previously unknown three-loop coefficient A^ 3 ' is obtained. The naive factorization of 

the cross section at small transverse momentum is broken by a collinear anomaly, which 
prevents a process-independent definition of xy-dependent parton distribution functions. 
A factorization theorem is derived for the product of two such functions, in which the 
dependence on the hard momentum transfer is separated out. The remainder factors 

^ into a product of two functions of longitudinal momentum variables and x\, whose 

renormalization-group evolution is derived and solved in closed form. The matching of 
these functions at small xt onto standard parton distributions is calculated at 0(a s ), 
while their anomalous dimensions are known to three loops. 



1 Introduction 

In collider processes with several disparate scales, fixed-order perturbative results become un- 
reliable since higher-order corrections are enhanced by large logarithms of scale ratios. The 
classic example of such a multi-scale process is the production of electroweak bosons with 
transverse momentum q? small compared to their mass M. The leading logarithmically- 
enhanced corrections in this kinematic region were resummed in [1, 2, 3]. An all-order formula 
for the cross section at small qr was obtained by Collins, Soper, and Sterman (CSS) [4], and 
explicit results for the ingredients necessary for resummation at next-to-next-to-leading loga- 
rithmic (NNLL) order were derived in [5, 6, 7]. The region of small qx is of phenomenological 
importance, since it has the largest cross section and is used to extract the W^-boson mass 
and width. In fact, the measurement of the charged- lepton qx spectrum now gives the most 
precise determination of Mw [8, 9, 10]. 

While the vector-boson q? spectrum is a classical example of an observable which exhibits 
logarithmic enhancements, analyzing its factorization properties is nevertheless rather subtle. 
In Section 2, we study the factorization properties of the Drell-Yan cross section at low trans- 
verse momentum in Soft-Collinear Effective Theory (SCET) [11, 12, 13]. At low qr, the cross 
section splits into two transverse-position dependent parton distribution functions (to which 
we will sometimes simply refer as "transverse PDFs") multiplying a hard function depend- 
ing only on the vector-boson mass M. However, it is well known that a naive definition of 
transverse-position dependent PDFs leads to inconsistencies due to rapidity divergences [14]. 
A proper definition requires introducing additional regulators beyond dimensional regulariza- 
tion. After taking the product of the two transverse PDFs these regulators can be removed, 
but finite terms depending on the hard momentum transfer q 2 = M 2 of the scattering process 
remain. We use the term "collinear anomaly" for this effect, since it describes a situation 
where a property of the classical theory cannot be maintained after including quantum cor- 
rections. We show that the anomalous terms have a very specific structure, which allows us 
to exponentiate the g 2 -dependent pieces in the product of two transverse PDFs. Once this is 
done factorization is restored. Specifically, we derive a factorization formula, which at fixed 
transverse displacement xt expresses the product of two transverse PDFs in terms of three 
functions capturing the dependences on the hard scale q 2 and the two light-cone momentum 
fractions £i and £2- The renormalization-group (RG) equations for these functions are de- 
rived and solved in closed form. This leads us to a consistent definition of the concept of 
XT-dependent PDFs, valid to all orders of perturbation theory. A simpler example for the 
occurrence of the collinear anomaly is the Sudakov form factor of a massive vector boson, 
which was discussed in SCET in [15]. The exponentiation of the associated InM 2 dependence 
was demonstrated in [16]. 

In the region qx 3> Aqcd, the transverse PDFs can be expanded in powers of x\. One 
then recovers the standard PDFs convoluted with xy-dependent kernel functions, which we 
compute at one-loop order in Section 3. A similar situation was studied in [17, 18]. Instead of 
the transverse PDFs relevant here, these papers considered a g + -dependent PDF, where q + is 
the small light-cone component of the partons emitted from the proton. What makes our case 
more complicated is the presence of the collinear anomaly, which neccessitates that both beam 
nucleons are considered at the same time and that an additional regulator is introduced. In 



our computations we use analytic regularization of the collinear propagators as in [19, 20]. We 
show that in xt space the In q 2 dependence induced by the collinear anomaly exponentiates to 
all orders in perturbation theory. In Section 4 we solve RG equations in SCET to derive, for 
the first time, a closed formula for the cross section in the region Aqcd ^ Qt ^ M directly 
in momentum (q^) space, which is free of large perturbative logarithms. This formula is the 
analog of the classical CSS formula, which is written in impact-parameter space, and it is valid 
to all orders of perturbation theory. 

Besides large logarithms, the perturbative series for the transverse-momentum distribution 
contains terms featuring a strong factorial growth at higher orders in a s , which must be 
summed to all orders to obtain a reliable result. In Section 5 we explain the origin of these 
terms and show how their summation can be implemented in closed form. In contrast to 
the CSS approach, our result for the resummed hard-scattering kernels is free of Landau- 
pole ambiguities, because we never perform scale setting inside integrals over the running 
coupling. Its connection to the CSS formula is discussed in Section 6. There we also derive 
all ingredients necessary for NNLL resummation, in particular also the previously unknown 
three-loop coefficient A^ 3 \ and show that it is not equal to the three- loop cusp anomalous 
dimension, as is often assumed in the literature (see e.g. [21, 22]). 

In previous analyses of qx resummation in SCET [23, 24, 25], which we review in Section 6, 
the collinear anomaly was regulated by keeping power-suppressed terms in the Lagrangian for 
the computation of the leading-power cross section. This complicates calculations, since then 
the various component functions no longer have homogeneous scaling in the SCET expansion 
parameter. More importantly, none of the previous works has addressed the resummation of 
the logarithms of q^/M 2 , which arise in the matching of the transverse-position dependent 
PDFs onto the standard PDFs. These define a class of next-to-leading logarithms, which 
cannot be resummed using the evolution of the hard matching coefficient. In the present work 
we show how these logarithms can be resummed to all orders in perturbation theory, working 
directly in momentum space using the formalism of [26]. Our work thus provides the first 
complete resummation of the large logarithms in the qx spectrum from effective field theory. 
A summary of our main results and conclusions are given in Section 7. Phenomenological 
applications of our approach will be presented elsewhere. 

2 Derivation of the factorization formula 

We study the Drell-Yan production of lepton pairs or electroweak gauge bosons with invariant 
mass M and transverse momentum q^ in the kinematical region where M 2 S> q\. Here q^_ 
denotes the transverse-momentum vector of the boson orthogonal to the beam axis, and we 
denote q\ = — q\ > 0. Our goal is to systematically resum large logarithms of the ratio M 2 jq\ 
to all orders in perturbation theory. Most of our analysis will assume that q\ ^> Aq CD is in 
the perturbative domain, but we will also discuss the case where q\ ~ Aq CD in some detail. 



2.1 Kinematical considerations 

For concreteness we consider the production of a lepton pair via a virtual photon with total 
momentum q^. The cross sections for Drell-Yan production of W and Z bosons can be obtained 
from the results presented here by means of simple substitutions summarized in Appendix A. 
We begin with the standard relation 

da = w~s W)* J d4xe ~ t9 ' x { ~ 9 ^ {Nl{p) N ' m rt{x) r{0) ]Nlip) N2{p)) ' (1) 

where 

J" = E Uw^q + gW^q) , 9l = g q R = e q (2) 

is the electromagnetic current. Keeping the left- and right-handed couplings separate has 
the advantage that our analysis can be carried over straightforwardly to the case of weak 
gauge-boson production. 

To analyze this process in SCET, we introduce two light-like reference vectors n and n 
satisfying n-n — 2, which are parallel to the directions of the colliding hadrons Aq and N 2 with 
momenta p and p. Any four momentum can then be split into light-cone and perpendicular 
components according to 

k» = n ■ k — + n ■ k — + k^ = kl + kt + k^ . (3) 

In the effective theory one defines a small expansion parameter A = qx/M and distinguishes 
fields whose momentum components (n-k, n-k, k±) scale differently with A. In the present case, 
the production of a lepton pair with transverse momentum q? requires one or more parton 
emissions into the final state, which balance that momentum. As a result, the substructure 
of the colliding hadrons is probed at distance scales of order 1/qr- The partons in the col- 
liding beam jets thus generically have so-called hard-collinear (he) or anti-hard-collinear (he) 
momenta scaling as 

p, c ~M(A 2 ,l,A), p^~M(l,A 2 ,A). (4) 

In SCET, one introduces different sets of hard-collinear and anti-hard-collinear quark and 
gluon fields describing the interactions of these partons [11, 12, 13]. While these fields have 
QCD-like interactions among themselves, two fields belonging to different sectors can interact 
only via the exchange of soft partons (only soft gluon interactions contribute at leading power 
in A), whose momenta scale like 

p s ~M(A 2 ,A 2 ,A 2 ). (5) 

Adding a soft momentum to a (anti-)hard-collinear momentum does not change its scaling 
properties. Since in our case the total transverse momentum of the hadronic final state must 
balance the transverse momentum q? = MX of the lepton pair, it follows that this final 
state must contain at least one (anti-) hard-collinear parton, and that soft fields give a power- 
suppressed contribution to its transverse momentum. As a result, we will see that the total 
contribution of an arbitrary number of soft emissions cancels in the final factorization formula 
for the Drell-Yan cross section at fixed q?. This is in contrast with the factorization formula 



for Drell-Yan production near threshold, in which soft modes give an important contribution 
[27, 28] (see [29] for a derivation using SCET). The absence of a soft contribution is an 
important fact. If there were a sensitivity to the soft scale fi s ~ q^/M, the Drell-Yan qT 
distribution would be non-perturbative even at high transverse momenta qx < a/^Aqcdj 
which would be a disaster for phenomenology. 

One might wonder whether also partons with momenta scaling as M(A, A, A) should be 
considered, which could contribute to the spectrum at leading power. On the one hand, the 
appearance of such a semi-hard mode would be quite unexpected from the point of view 
of SCET, since it is not allowed to interact with the (anti-)hard-collinear partons. Indeed, a 
semi-hard mode does not contribute to any other collider-physics process analyzed with SCET 
so far. One the other hand, because the transverse momentum of the final-state partons is 
restricted to scale like MX, it is not obvious that the emission of semi-hard partons is really 
irrelevant. In Section 3, we will show that the semi-hard region does not contribute if the 
collinear anomaly is regularized analytically. In this respect our analysis differs from a recent 
study in [25]. We will comment on this paper in more detail in Section 6. 

In SCET, the current (2) is matched onto (we adopt a regularization scheme with anti- 
commuting 75) [30, 31] 



J^ ^C v (-q 2 -te^^J^^ix-cSti^-^S^kc + glx-cStr^^S^hc) 



(6) 



where the matching coefficient Cy depends on the hard momentum transfer Q 2 = —q 2 = —M 2 . 
For time-like processes this coefficient is complex. From now on we will suppress the — is 
prescription that defines the sign of the imaginary part of Cy. The effective fields Xhc = Wl c 61c 
and Xhc = Wj-£fa are the usual gauge-invariant combinations of effective (anti-)hard-collinear 
quark fields and Wilson lines. 1 They satisfy fiXhc = and fix~hc~ = 0- These fields are obtained 
after the SCET decoupling transformation has been applied, which removes the interactions 
between soft and (anti-)hard-collinear fields in the leading-order SCET Lagrangian [12]. This 
introduces the soft Wilson lines S n and S n into the effective current shown above. In [29], we 
have derived a factorization theorem for the Drell-Yan process in the threshold region. Many 
more details about the effective-theory formalism can be found in this reference. 

Using a Fierz transformation along with some elementary Dirac algebra, and averaging 
over nucleon spins, the hadronic matrix element in (1) can be rewritten as 

\9l\ 2 + \9 q R ? 



-g, v ) (N 1 (p)N 2 (p)\ Jrt(x) r(0) \N 1 (p)N 2 (p)) -► \C v (-q 2 ,fi)\ 



2 



^ 2N, 



c 



x Wv Y (x) (Aq(p)| Xhc (x) \ XhM \NM) (N2(P)\ XttM \ Xhix) \N 2 (p)) 



(7) 



where 

Wv Y (x) = -^(0\Tr\T(Si(x)S n (x))T(Sl(0)S n (0))]\0) (8) 



N, 



c 



is a soft Wilson-line correlator. Up to this point our discussion is completely analogous to that 
in [29]. We now perform the multipole expansion appropriate for the kinematical situation 



1 As usual in SCET we restrict ourselves to gauge transformations that vanish at infinity. 



considered here, i.e., we expand the fields in derivatives corresponding to suppressed momen- 
tum components [13, 32]. Here we encounter a difference with our previous analysis, which is 
due to the different kinematical requirements on the transverse momentum of the photon. In 
the present case the separation of the fields scales as x ~ M _1 (l, 1, A -1 ), which is conjugate to 
the hard photon momentum q ~ M(l, 1, A). Since derivatives on the soft fields scale as A 2 , the 
leading term in the expansion is obtained by evaluating the soft Wilson lines at x — 0. In the 
(anti-)hard-collinear fields, the dependence on x± and on the light-cone components conjugate 
to the large momentum components is unsuppressed and must be kept. The expanded result 
then takes the form 

x (JVi(p)| Xhc(x+ + zi.) I X*c(0) |iVi(p)) (N 2 (p)\ x E (0) I xc(i- + zi) |JV 2 (j>)> . 

The definition (8) implies that Wdy(0) = 1, and hence the soft contribution cancels out in the 
result for the cross section. Physically, the reason is simply that soft gluons carry transverse 
momenta of order q^/M, which for M ^> qx implies that an arbitrary number of soft gluons 
can be emitted into the final state without changing the total transverse momentum q?. Hence, 
real and virtual soft-gluon effects cancel by the KLN theorem. 

We recall at this point the definitions of the parton distribution functions (PDFs) [33, 34, 
35]. In terms of SCET operators, they read 

(f> q/N (z,u) = ±- I ' dte~ izt ^ (N(p)\ xitn) | x (0) \N(p)) , 

(10) 

Zn ' V I j + -iztn-p 



<f> g/N (z,n) = ^fj dte-* ztn - p (N(p)\ -A^(tn)A>i(0) \N(p)} , 

where < z < 1, and for simplicity we suppress the he labels on the fields. The nucleon iV 
carries momentum p along the n direction. The anti-quark PDF <f)q/^(z,fi) is given by the 
same matrix element as in the first line, but with the sign in the exponent reversed. Because of 
the appearance of the transverse displacement vector x±, the hadronic matrix elements in (9) 
are not of the form of the matrix elements of light-ray operators defining the PDFs. Instead, 
one defines the generalized, xy-dependent PDFs (with x\ = — x\ > 0) [34, 35] 

B q/N (z, 4, p) = ±Jdt e-^ (N(p) | xitn + x ± ) | x(0) \N(p)) , (11) 

and similarly for the gluon and anti-quark cases. Their Fourier transforms with respect to 
xt are referred to as transverse-momentum dependent PDFs. Naively, then, the differential 
cross section (9) can be expressed in terms of a convolution of the hard matching coefficient 
|CV(-M 2 ,/x)| 2 with x T -dependent PDFs, 

d 3 a Alia 2 >~, , r0 s 1 2 1 f , 9 ,-„ _ 

\C V (-M 2 ,n)\ — d 2 x ± e- %q± - x± 



dM 2 dq 2 dy 3N c M 2 s ' VK '^ y| 4vr 



E 



x > 4 



Bq/m (£1, x 2 t, (*) &q/N 2 (6, x 2 T ,fi) + (q^ q) 



\M 2 



(12) 



where we have used that g\ — g q R — e q , and have defined 

6 = v^, 6 = v^e^, with T = < = M1±A, (13) 

s s 

We have also employed that d A q 9(q°) 5(q 2 —M 2 ) = | d 2 q± dy = | dg T dy, where the last identity 
holds after integration over the azimuthal angle. The above formula appears to achieve the 
desired factorization of the hard and hard-collinear scales, M 2 and q\ ~ x^ 2 . 

2.2 Collinear anomaly and refactorization 

The formal derivation of factorization just presented is spoiled by quantum effects. This can 
be seen from the fact that the RG equation for the hard matching coefficient contains a term 
proportional to a "cusp logarithm" of the hard scale q 2 = M 2 [26], 



C v (-q 2 ,p) 



A 4 



CV(-g 2 ,/i). (14) 



dlnfi 

Its coefficient r^ s p (a s ) is the cusp anomalous dimension in the fundamental representation. 
The quantity / ~f q (equal to 7 y /2 in [26]) refers to the quark anomalous dimension as defined 
in [36, 37]. Here and below, the coupling a s without an explicit scale argument always refers 
to a s (fi). RG invariance of the physical cross section (9) requires that the evolution equation 
for the product of the two transverse PDFs must contain the same ln(q 2 /fi 2 ) term as in (14), 
but with the opposite sign. This fact is incompatible with the concept of universal (i.e., 
process-independent) transverse-position dependent PDFs. As we will see in Section 3.1, the 
XT-dependent PDFs as given in (11), in which gauge invariance is ensured by means of light- 
like Wilson lines, are not well-defined in dimensional regularization and require an additional 
regularization of light-cone singularities. The anomalous q 2 dependence arises because the 
formal factorization of hard-collinear and anti-hard-collinear fields in the SCET Lagrangian 
in the absence of soft interactions is invalidated as soon as one introduces a regulator to give 
meaning to the loop integrals appearing in explicit evaluations of SCET diagrams. We will 
refer to this effect as the "collinear anomaly". Only the product of two transverse PDFs 
referring to hadrons moving in different directions is regularization independent. A similar 
phenomenon has been encountered in the context of SCET analyses of the B — y it form factor 
[38, 39, 40], where however the factorization breakdown happens at non-perturbative scales. 2 
The collinear anomaly is, of course, not a new quantum anomaly of QCD. Rather, it is 
a feature affecting certain matrix elements in SCET - the effective theory of QCD relevant 
to the derivation of QCD factorization theorems. In SCET this effect is an anomaly in the 
usual sense that quantum corrections destroy a symmetry of the classical theory. At the 
classical level, the hard-collinear sector in the effective Lagrangian does not know about the 
anti-hard-collinear momentum p, and it is thus invariant under the rescaling transformation 
p —y Xp. Likewise, the anti-hard-collinear sector in the effective Lagrangian is invariant under 
the rescaling transformation p —y Xp. We will show in Section 3.1 that in the presence 



2 Aftcr submission of this paper, we learned that in this context the term "factorization anomaly" was used 
in lectures delivered by M. Beneke (http://theor.jinr.ru/ hq2005/Lectures/Beneke/Beneke-Dubna-05.pdf). 

6 



of the regulators required to give meaning to SCET loops integrals only the subgroup of 
transformations with AA = 1 is left unbroken. While q 2 = 2p ■ p is not invariant under 
independent rescalings of p and p, it is invariant under the subgroup A A = 1. This explains 
why an anomalous dependence on q 2 can arise in the regularized theory. Fortunately, the 
anomalous terms have a very specific structure, and this will allow us to exponentiate the q 2 
dependence of the product of two transverse PDFs, thereby restoring factorization. 

Because of the dependence of the transverse PDFs on the hard scale of the underlying 
process, relation (12) does not accomplish a complete separation of the hard and hard-collinear 
scales q 2 = —M 2 and q\. It is thus not yet a useful factorization formula. In order to complete 
the factorization and resum all associated large logarithms, it is necessary to control the q 2 
dependence of the product B q /^ Bq/N 2 °f transverse PDFs to all orders in perturbation theory. 
We will show in Section 3.2 that in xt space this product can be refactorized in the form 

/ x 2 2 \ -Fqqi^T^) 

[B q/Nl (z u x 2 T , id) Bq /N2 (22, x 2 T , //)] q2 = I 4e ^ 2/yg J B q/Nl ( Zl , x 2 T , fl) B q/N2 (z 2 , x 2 T , fl) , 

(15) 
where the exponent F qq depends only on the transverse coordinate and the renormalization 
scale. The bracket on the left-hand side of (15) indicates the hidden q 2 dependence induced by 
the collinear anomaly. The functions -Bj/at on the right-hand side are independent of the hard 
momentum transfer. While we do not have an operator definition of these functions, they are 
uniquely defined by relation (15). All q 2 dependence is now explicit and controlled by the 
function F qq . Note that for \x ~ x^ 1 the g 2 -dependent prefactor resums all large logarithms of 
the hard scale, while F qq (xj^,fj.) has a perturbative expansion in a s (fi) with 0(1) coefficients. 
RG invariance of the cross section (9) implies the evolution equations 



dr qq {x T) fi) _ oinF 
d 



dkifi - 2 W«s) 



dhifj, 



B q/N (z,x T ,fi) 



rf usp K)in^-2 7 v, 



(16) 



B q / N (z,a%.,(i) 



The first relation completely determines the scale dependence of F qq . Note that the structure 
of the evolution equation for Btm is completely analogous to that for the hard matching 
coefficient in (14). As a side remark, let us add that for gluon-initiated processes such as 
Higgs-boson production analogous evolution equations hold for the quantities F gg and B g / N , in 
which r^ sp is replaced with the cusp anomalous dimension T^ in the adjoint representation, 
and 7 9 is replaced with the gluon anomalous dimension / ~f 9 as defined in [36, 37]. Formal 
all-order solutions to the evolution equations (14) and (16) will be presented in Section 4 and 
Appendix B. Moreover, we will conjecture in Section 3.2 that the exponent F qq in (15) is 
constrained by the non-abelian exponentiation theorem [41, 42], which implies the Casimir- 
scaling relation 



up 0/ 



at least to three-loop order. 



From the point of view of our SCET analysis, equations (15)— (17) contain all there is to 
say about the concept of XT-dependent PDFs. The factorization theorem (15) specifies in a 
precise, all-order way that the dependences on the variables q 2 , zi, and z 2 can be factorized at 
fixed x\ into three functions. The exact evolution equations for these functions are given in 
(16) and can be solved in closed form. Finally, relation (17) relates the exponents of the in- 
dependent terms in the quark- ant iquark and gluon-gluon channels. We will show in the next 
section that for small transverse separation, x\ <C Aqq D , the exponent F qq - and the functions 
-Bj/at can be evaluated in perturbation theory, the latter ones being related to standard PDFs 
in a calculable way. For larger values x\ ~ Aq^ d , however, the functions F qq and -Bj/at are 
genuine, non-perturbative objects, which must be modeled or extracted from fits to data. 
Relation (15) provides a rigorous basis for such modeling. We will discuss below how an 
experimental determination of the non-perturbative functions could be performed in practice. 

A review of the conventional view on transverse-position dependent PDFs can be found in 
[14] (see also [43]). There it is emphasized that these objects are indeed problematic, because 
they are subject to light-cone singularities not accounted for by standard renormalization 
methods. These singularities originate from gluons with very large rapidity relative to the 
colliding hadrons, and regularizing them requires some sort of a rapidity cut-off. While in the 
present work we use analytic regularization for this purpose, two different proposals were made 
in [14]: the introduction of Wilson lines off the light cone (and subject to certain constraints) 
[45] , and a "generalized renormalization prescription" consisting of multiplying the original xt- 
dependent PDFs (11) with a gauge-invariant factor that cancels the extra divergences [46, 47]. 
These constructions have been explored at one-loop order only. Most practical applications 
of xy-dependent PDFs employ a definition based on an axial gauge v ■ A — 0, as originally 
introduced by Collins and Soper in [34, 35] and denoted by Vi/N(z,XT,lJ>',C)- This introduces 
a dependence on a "gauge parameter" ( = (2p ■ v) 2 /v 2 and thereby a sensitivity to the energy 
p° of the external nucleon in the reference frame defined by v° = l. 3 The product of the 
gauge parameters associated with the two nuclei obeys v^CiG = 2p ■ p = a 2 /(ziz 2 ), which 
is a gauge-independent result [47]. In this way the hard momentum transfer q 2 enters in 
the Collins-Soper approach, even though the precise way in which the ( dependence cancels 
between the XT-dependent PDFs and their hard function is not very transparent. The reason 
is that the ( dependence of Pj/at(z,xt, fJ>; C)> as given by relation (6.23) in [34], is much more 
complicated that the q 2 dependence in our result (15). 

Despite of these efforts, the precise concept of transverse-position dependent PDFs has 
remained elusive [49]. The insight that only the product of transverse PDFs relevant in a 
given process is physical and can be defined in a regularization-independent way is a crucial 
new ingredient of our analysis. Moreover, our results (15) and (16) define the functions B^m 
and F q q and their evolution equations even at a non-perturbative level. On the other hand, 
it is an interesting open question whether the functions Bi/n are universal and appear also 
in other processes (see [50] for a related discussion). It seems conceivable to us that in more 
complicated processes involving more than two directions of large energy flow in the hard- 
scattering event, relation (15) could take a more complicated form. 



3 In the equivalent description using Wilson lines pointing in a time-like direction v, the gauge dependence 
is replaced by a dependence on the auxiliary vector v (see e.g. [48]). 



In terms of the objects defined in (15), the differential Drell-Yan cross section at fixed 
invariant mass M, transverse momentum q?, and rapidity y of the Drell-Yan pair can be 
expressed as 



d 3 a 


Ana 2 


\C ( M 2 u)\ 2 l fd 2 x :-*"■*' (^^ 


\ -F qi (x%,n) 


dM 2 dq\ dy 


3N c M 2 s 

x E e ? 

q 


\c V { m ,fi)\ ^ j dx ±c y Ae -2 1E 

B q / Nl (£i, Xt, h) Bq /N2 (&, Xt, h) + (q <H> q) 


J 

\M 2 



The disparate scales M 2 and q\ ~ x^, 2 are now completely separated, and all large logarithms 
can be resummed by choosing /j, ~ q T (or jj, equal to a few GeV in the case where qx ~ Aqcd) 
and employing the RG-improved expression for CV(— M 2 ,/z) given in relation (52) below. 
Corrections to the leading term in the factorization formula are suppressed by powers of the 
ratio q\jM 2 <C 1. Also, as written above, the formula holds irrespective of whether or not the 
transverse momentum is a perturbative scale. Taking a Fourier transform of the cross section, 
it is possible to get direct access to the x-r-dependent PDFs as given in the factorization 
theorem (15). We find 



9M 2 s 
Ana 2 



dqr Jq(QtXt) 



d 3 



a 



o 



dM 2 dql dy 



\C V (-M 2 ,fi)\ 



x 2 T M 2 

4 e -2 7 B 



-Fqq( X T>V) 



X 



B q /N 1 {^l-,X T ,H)B 1 



y q/N 2 (C2,X T ,IJ,) 



(q++q) 



o 



x 2 T M 2 



(19) 



q 

By varying Xt, M 2 , s, y, and the beam nuclei JVi, N 2 , one can (at least in principle) map 
out the functional dependences of F qq and certain combinations of transverse PDFs on x\ and 
£i, much in the same way as the standard PDFs are constrained from fits to Drell-Yan cross 
sections. While for xt «C Aqq D the right-hand side of (19) can be calculated in terms of known 
PDFs (see below), for xt ~ Aq^ d the above relation provides access to the non-perturbative 
behavior of F qq and of the transverse PDFs. This can help to constrain phenomenological 
models of these functions, which are needed e.g. for a precision determination of the mass 
of the W boson. We emphasize that the above relation only holds for x\ ^> 1/M 2 , because 
otherwise the power corrections to our factorization formula become large. It can therefore 
not be used to study the xt — > limit of the functions F qq or Bi/ N . 



2.3 Simplifications at large q^ 

For given transverse momentum q?, the Fourier integral in (18) receives important contri- 
butions from transverse separations Xt < q? 1 only For large transverse momenta in the 
perturbative domain, q\ 3> Aq CD , we therefore need the xy-dependent PDFs at transverse 
separation xt <C A Q q D . In this case these functions obey an operator-product expansion of 
the form [4, 34, 35] 

11 dz . 



Bi/n(£,,x 



T,A*) 



£ 



l *i-3 



>, x T , fj) (j) j/N (£/z, //) + C(A 



2 
QCD 



X T ) 



(20) 



In the context of SCET, generalized PDFs defined in terms of hadron matrix elements in 
which collinear fields are separated by distances that are not light-like are referred to as beam 
functions. For such functions an analogous expansion was considered in [17], and an expression 
for the one-loop kernel of the quark beam function was derived in [18]. The evolution equations 
for the new kernels Li^j follow when we combine (16) with the standard DGLAP equations 



d 



din jj, 



/■i i 

■ J z u 



(21) 



We obtain 



din// 



l q<r- 1 



Z, x T , n) 



X Li 

r c F uspK)ln^^ 



2 7 > s ) 



L q<—t 



z, x T , n) 



r 1 d 
~^2 —Ig^ j (u,x^,fi)V j ^ i (z/u,fi). 

A J Z U 



(22) 



Because of the complicated form of the DGLAP equations, no closed solution can be derived. 
Neglecting power corrections of order A.Q CD /q^, we can use relation (20) to express the dif- 
ferential cross section (18) as a convolution of perturbative, factorized hard-scattering kernels 



9 1 

2 S/f2 ,a \r~i I t,j2 ,.\l z ■ / j2_ „-«<J_L-a:_L 



C qq ^ ij (z 1 ,z 2 ,q T ,M ,/j,) = \C V (-M ,ti)\*-^- / ^ Xl - e 
with ordinary PDFs. The result reads 



x 2 T M 2 ^ ~ f ^ x t^ 
4 e -2-y E 



(23) 



d 3 



a 



47TO 2 



dM 2 dq\ dy 3N c M 2 s 



E4EE 



x 



a 



i=q,9 3=<i,g 

.2 n/rl 



1 dz\ f 1 dz 



fc Zi J& z 2 



(24) 



qq->v 



z 1} z 2 ,q T ,M ,/j)(f>i/Ni^i/zi,tJ)<f>j/N2^2/z2,fi) + (q,i*+ q,j) 



This formula, as well as relations (25) and (27) below, receive power corrections in the two 
small ratios q\jM 2 and AQ CD /g|,. This will not be indicated explicitly. 
Integrating this result over rapidity, with |y| < ln(l/r), we obtain 



d 2 a 



Ana" 



9:yj y>9 z\Z2> 



dz\ dz 2 

Z\ z 2 



dM 2 dq\ 3N c M 2 s 

x C gg _>ij(zi, z 2 , q 2 T , M 2 , //) ff i:j ( , /i) +(g,i« q,j) 

\Z\Z 2 / 



(25) 



where the parton luminosities are defined as 

1,1 dz 



ffij(u,V>) 



(f) i/Nl (z,fj,)(f) j/N2 (u/z,n) 



(26) 
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Figure 1: Matching of an analytically- regularized QCD graph onto SCET diagrams. 



It will also be useful to study the total cross section defined with a cut qx < Qt, which vetoes 
single jet emission. Neglecting the dependence of the variable r in (13) on gf., which is a 
power-suppressed effect, we obtain from (25) 

1=1,9 0=1,9 



da 



dM 2 



„ <n 3N c M 2 s 

qrSQr q 



z 1 z 2 >M 2 /s 



X 



min(Qy , z±Z2S—M 2 ) 



dq 2 T Cgg^Zi, z 2 , q 2 T , M 2 , fj) ff tj , /j,) + (q,i*+q, j) 

\Z1Z2S / 



3 Calculation of the kernels I q ^ q and I q ^ g 

We now perform a perturbative calculation of the relevant kernels I^j entering the factor- 
ization formula (23) at first non-trivial order in a s . Since we do not have explicit operator 
definitions of the refactorized transverse distribution functions B^m, we analyze instead the 
original functions Bim defined in (11), keeping in mind that only products of two such func- 
tions referring to different hadrons are well defined. If we write an operator-product expansion 
analogous to (20) 

Bi/ N (^,XT,n) = ^2 — z i^j( z ^ x T^) ( t ) j/N^/z,fi) + 0(A 2 QCD x 2 T ), (28) 

j J t z 

it follows that the products of two Zj«-j functions are well defined and obey a factorization 
formula analogous to (15). 

3.1 One-loop results 

Perturbative expansions for the kernels 1-i^j can be derived from a matching calculation, in 
which the matrix elements in (10) and (11) are evaluated using external parton states carrying 
a fixed fraction of the nucleon momentum p. The tree-level result is obviously given by 



Ia-j(z, Xt, (j) = S(l - z) Sij + 0(a s ) . 



(29) 



However, when trying to evaluate the one-loop corrections, one finds that they are ill-defined 
in dimensional regularization due to light-cone singularities. To give meaning to the corre- 
sponding loop integrals requires introducing additional regulators. The simplest possibility 
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is to employ analytic regularization, as is common in the context of asymptotic expansions 
[19, 20]. In the context of SCET this method has been used in [15, 39]. One starts by re- 
considering the QCD diagrams contributing to the process and raises all propagators through 
which the external hard-collinear momentum p flows to a fractional power, 



,2a 



_(p _ k )2 _ l£ ~* [_ (p _ k) 2 _ i£ ]l+a > ( 30 ) 

and similarly for the anti-hard-collinear propagators, but with a different regulator f3 and an 
associated scale v-i- For QCD diagrams, such as the first graph in Figure 1, the modification 
is trivial in the sense that the limits a — >■ and /3 — > are smooth as long as the dimensional 
regulator d = 4 — 2e is kept in place. However, with analytic regulators the contributions of 
the different momentum regions are now well-defined individually and one can check which 
regions give non- vanishing contributions to the expansion of the loop integrals. One finds that 
only the hard, hard-collinear, and anti-hard-collinear regions contribute. If the diagrams are 
evaluated off-shell, then also a soft contribution arises in the individual diagrams, but one 
easily verifies that the contribution of the semi-hard mode vanishes. Since such contributions 
were considered in the literature, we now explicitly show that they are absent in analytical 
regularization. To do so, one assumes that the gluon momentum k in the QCD diagram in 
Figure 1 scales as M(A, A, A) and then expands the diagram in powers of A. At leading power, 
the relevant phase-space integral becomes 

d d k j— T - B 5{k 2 ) 9(k°) e v-*-*±-*± . (31) 

(n • k — ie) (n • k — ie) 

After the expansion, the integrand involves the usual eikonal propagators characteristic for 
soft emissions, which are raised to fractional powers because of the analytic regulators. It 
is important that not only the propagator denominators, but also the Fourier exponent is 
multipole expanded: 

(p - k) ■ x = p ■ x - k± ■ x ± + 0(A) , (32) 

which follows from the scaling x ~ (1, 1, A -1 ) derived in Section 2. Performing the integration 
over the light-cone components of the gluon momentum, one finds that the integral (31) is 
scaleless and vanishes. The same argument applies to multi-loop integrals involving semi-hard 
propagators. 

We note that the multipole expansion was not performed in [25], which explains why similar 
integrals were found to be non-vanishing in this reference. However, the expansion is a crucial 
ingredient to achieve scale separation for effective theories in dimensional regularization. It 
is equally important in the strategy of region technique [19, 20]. Without performing the 
expansion integrals pick up contributions from several regions, and care needs to be taken to 
avoid double counting. 

Having shown that the soft and semi-hard regions do not contribute, let us now turn to the 
hard-collinear contributions. Since the original diagrams are well defined without analytical 
regulators and are obtained by adding up the contributions from the different regions, we are 
guaranteed that the limits a — > and /3 — > can be taken in the sum of all diagrams and 
that the final result is independent of the regularization scheme. Individually, however, the 
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Figure 2: One-loop diagrams contributing to the matching coefficients I qi -q (top row) and 
Iq^g (bottom row). The vertical lines indicate cut propagators. 

diagrams in each sector involve divergences in the analytical regulators. If the momentum 
k in (30) is hard-collinear, as in the first SCET diagram in Figure 1, the a dependence in 
the effective theory takes the same form as in QCD. If, on the other hand, the momentum 
k is anti-hard-collinear, then the propagator is far off-shell and in SCET is represented by a 
Wilson line, as shown in the second diagram in Figure 1. Using the replacement rule (30) and 
performing the appropriate expansions, we find that the Feynman rule for a gluon emission 
from the anti-hard-collinear Wilson line Wj^ in the current operator (6) gets replaced by 



n" 



n ■ k 



->■ 



v\ a n^n ■ p 



ie 



[n 



kn ■ p 



IE) 



,1 + Q 



(33) 



Note that, as mentioned earlier in Section 2.2, the regularized Feynman rule for the anti-hard- 
collinear Wilson line is no longer invariant under the rescaling transformation p — > Xp. As 
seen in Figure 1, the regulator a plays a double role: it regularizes the fermion propagators 
in hard-collinear diagrams and the Wilson lines in anti-hard-collinear diagrams. Both classes 
of diagrams develop singularities in the limit — > followed by a — > or vice versa, which 
cancel in the sum of the results from both sectors. 

With the regularization in place, let us now turn to the evaluation of the one-loop correc- 
tions to the kernels X q ^ q = Tq-^q. The relevant effective theory diagrams are shown in the first 
row of Figure 2. There is no need to consider diagrams with external-leg corrections on only 
one side of the cut, because these give identical contributions to £>j/jv and 4>i/N and thus do not 
change the tree- level result (29). Working in Feynman gauge, we find that the contribution of 
the first diagram is finite in the limit where the analytic regulators are sent to zero as long as 
e is finite. A divergence arises only in the expansion in e, after which we obtain 



^■q^q\ Z i X Ti^) 



C F a s 
2vr 



'1 



1 



+ Lj 



In 



OC'-pLl 

4 e -l7E 



(34) 



up to (9(e) terms. As before, a s = a s (/i) always refers to the running coupling evaluated at 
the scale /i, unless indicated otherwise. Moreover, /i denotes the renormalization scale defined 



in the MS scheme. The fourth diagram gives a vanishing result, X\ 



q^q 



0, but the remaining 
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two diagrams in Figure 2 are non-zero and only well-defined with analytic regulators. Both 
give the same result, and for their sum we obtain 

a*.*,) = %v™ kv ay ,, 2 > T £:°} (*pr. (35) 



j-b+c 



2tt 



W/ W/ (i-^) 1- "^ r(i + «) V 4 ; 

Like in full QCD, the analytic regulators must be taken to zero before taking the limit e — >■ 0. 
The result depends on the order in which the limits a — > and /3 — > are performed. Ex- 
panding first in /3 and then in a, the light-cone singularities are regulated by the a parameter, 
and we find for the sum of all four one-loop diagrams 



Zj4~q(Z, x T , n) 



C F a s 

2tt 



+ L_ 



+ 8(l-z) 



2 

a 

2 



2 In 



/// 



8(1 -z) + 



1 + z' 



+ L 2 + 



7T~ 

If 



'1-z) 



(36) 



If the expansions are performed in the opposite order, then /3 acts as the analytic regulator, 
and we obtain 



z 



q<—q\ Z i X Ti A 4 / 



P reg. 



CfCXs 

2tt 



1 



+ L A 



2 q 2 

-- + 21n^ 

P vi 



8(1 -z) + 



1 + z 2 



(37) 
The above results refer to the kernel associated with hard-collinear partons, which propa- 
gate along the n direction. Let us now consider what happens when we calculate the corre- 
sponding kernel for anti-hard-collinear fields. In this case we get the same answer but with 
a, v\ and /3, v 2 interchanged. We then find that in the product of a hard-collinear and an 
anti-hard-collinear kernel function the analytic regulators disappear, no matter in which order 
the limits a — > and /3 — > are taken. This product is thus regulator independent and well 
defined in dimensional regularization. After MS subtractions, we obtain 

\Zqi-q(Zl, X T , fijZq^q{Z2, X T , fJ,)\ g2 

Cf&s 



S(l- Zl )S(l-z 2 ) 
C F a s 



1- 



2tt 



7T 

~6~ 



2ty 



8(1 - Zl ) 



i_ + 4 

1-32 



2L ± ln^- + L 2 ± - 3L ± + 

H 2 

- (1 - z 2 ) + ( Zl ++ z 2 ) 



(31 



0(al 



Next we calculate the kernel X g4 _ 9 at one-loop order. This function vanishes at tree level, 
and at one-loop order it follows from the evaluation of the diagram shown in the second row of 
Figure 2. There is no need for analytic regularization in this case, and after MS subtractions 
we find 

Z q ^ g (z, 4, /i) = -^ {L ± [z 2 + (1 - z) 2 ] - 2z(l -z)} + 0(a 2 ) . (39) 

From (38) and (39) we can extract the one-loop expressions for the renormalized function 
F q q and the renormalized kernels /j<_j relevant for Drell-Yan production. We obtain 

C F a s 



F q q(L ± ,a s ) 



L 



TV 



0(a 



2 ) 

Sl 1 



(40) 
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and 



I qi - q (z,L ± ,a s ) = 5(1 - z) 



C F a s 



1 + 



C F a s 
Aix 



7X 



2\ 1 



L\ + 3L ± - — 



L Q^g 



[z,L ± ,a s 



2n 
T F a & 

2vr 



£± W*) - (1 



L±P q <- g (z) - 2z(l - z) 



6 , 

- c(« s 2 ) , 



(41) 



with Lj_ as denned in (34). The kernel I q ^ q is given by the same expression as I q -+ q , and 
I q -+ g has the same form as I q ^ g - Note that, with a slight abuse of notation, we have changed 
the arguments x\ and ft in F qq and the kernel functions to L±_ and a s , as this will be more 
convenient from now on. In the above expressions 



*W*) 



1 + 2 2 
1-2 



w*) 



2 2 + (l 



(42) 



are the one-loop Altarelli-Parisi splitting functions, defined as 



' qi-q\ x i A* J 



C F a s 



TX 



P q ^Jx) + 0(a 2 s ), V q ^ g (x,v) 



T F a s 



7X 



P q ^ g (x) + 0(a 2 s ) . (43) 



It is straightforward to check that our one-loop results (38) and (39) satisfy the general evo- 
lution equations (22). 

3.2 All-order dependence on the hard momentum transfer 

The appearance of a logarithm of the large momentum transfer q 2 in the matching condition 
(38) appears strange at first sight, since it arises from the evaluation of hard-collinear and anti- 
hard-collinear loop graphs in the effective theory, in which these two sectors are decoupled 
from each other at the Lagrangian level. Naively, we would thus expect a dependence on 
the scale x F only. For \x of order a typical hard-collinear scale the resulting logarithm is 
parametrically large, so that a s \n(q 2 /fi 2 ) ~ 1 in RG power counting. The question then arises 
if higher powers of such logarithms appear in higher orders of perturbation theory, and if this 
is the case, how these logarithms can be resummed to all orders. 

A simpler example for the occurrence of the collinear anomaly is the Sudakov form factor of 
a massive vector boson, which in SCET has been discussed in [15, 16]. The similarity with our 
case arises because the (^-function requiring that the square of the total transverse momentum 
k±tat of all cut propagators in the discontinuity of the forward-scattering amplitude be equal 



4) 



which indeed 



to q\ can be rewritten as the discontinuity of a "propagator" l/(fcj_tot 
looks similar to a massive propagator. Using arguments along the lines described in [16], we 
will now argue that the logarithm of the product of two transverse PDFs B.i/n 1 Bj/n 2 , or of 
two kernels Xi^k^ji-i, is linear in ln(q 2 /fi 2 ) with a universal coefficient —Fij(x^, //) that is 
independent of z±, Z2 and k, I, as shown in (15). This result implies that the large logarithms 
exponentiate in x F space. 

Because of the formal decoupling of the hard-collinear and anti-hard-collinear sectors in the 
SCET Lagrangian, the product Bi/^ Bj/n 2 is a product of the sum of all hard-collinear graphs 
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times the sum of all anti-hard-collinear graphs, each defined by means of analytic regulators. 
Let us denote the logarithms of the sums of all (anti-)hard-collinear graphs by 



/ u 2 q 2 \ 

\nB i/Nl = ln(J^ he graphs) = %(fi, L±,a a (/j,), In — , In — 

In Bj/ N2 = ln( V" he graphs) = Hj ( £2, L±, a s (/i), In — , In — J 



(44) 



The need for analytic regularization introduces logarithmic dependence on the auxiliary scales 
V\ and z/ 2 associated with the a and /3 regulators, see Section 3.1. The precise form of the 
results for Hi and Hj (but not for their sum) depends on the order in which the limits a — > 
and /3 — > are taken. All that matters is that the same prescription is used throughout the 
calculation. In the hard-collinear sector the logarithmic dependence on v\ can for dimensional 
reasons only appear in the form of ln(/z 2 /z/ 2 ), since no other momentum scale is introduced by 
the a regulator. In the anti-hard-collinear sector, on the other hand, the induced logarithmic 
dependence on v\ is accompanied by n ■ p. By Lorentz invariance, a logarithm of the "foreign" 
momentum component n-p can only appear in the form of n-pn-p = q 2 , so that we encounter 
ln(g 2 jv 2 \ For the /3 regulator with associated scale z/2, the role of the two sectors is reversed. 
This explains the dependencies on the two regulator scales shown in (44). 

The requirement that the product of the two transverse PDFs must be independent of the 
two regulator scales implies that 

— ^- {Hi + Hj) = = -^— {H t + H J ) . (45) 

These conditions enforce that Hi and Tij are linear in their last two arguments with coefficients 
that are independent of £i and £ 2 , and hence we are free to write 

lnB i/Nl B j/N2 = Hifa, L±, a a (fi)) + Hj{&, L±, a s {fi)) - Fij(L±, a fl (/x)) (in — + L± J . (46) 

We have used the freedom that the decomposition is unique up to a function of L± to make the 
coefficient of F q q scale independent. With the identification Hi = \nBi/^ 1 and Hj = In Bj/n 2 
this proves relation (15), where we have considered the special case i — q, j — q. 

Let us now collect what can be said about the function F q q based on general principles. 
Generalizing our one-loop result to higher orders, we can write the perturbative expansion of 
F q q in the form 

oo 

Fqq(L u a s ) = J2dUL^(^y, (47) 

71=1 

where d\(L±) = 4CfL±. The first evolution equation in (16) then implies the recursion 
relation 

n— 1 

dl'(L ± ) = 1^ + J2 mPn-i- m dl(L ± ) , n > 1 , (48) 

m=l 
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where the prime denotes a derivative with respect to L±, and as usual we have expanded the 
cusp anomalous dimension and f3(a s ) = fida s /dfj, as 



\Oi s ) 



E r »-.(£) ■ ««.) = -2<».Ea-'(S) ■ (49) 

n=l n=l 

For the first two expansion coefficients, we obtain 



d\(L ± ) = r*L ± + d\ , d»(L ± ) = V J^Ll + TfL ± + d\ , (50) 

where d q n = rf^(0) with g^ = 0. The expansion of the corresponding function F gg for Higgs 
production can be written as in (47) but with coefficients d g nl which obey analogous equations 
in which r^ sp is replaced by r^ sp . We will later discuss how the two-loop coefficients d\ 9 
can be extracted from existing calculations of higher-order corrections to Drell-Yan and Higgs 
production cross sections derived in fixed-order perturbation theory [6, 7]. The result is 

d\ d% „ /808 _"\ 224 



trt =c "{iT- 28C: ')-w T ^- (51) 

These coefficients contain only maximally non-abelian color structures. This leads us to con- 
jecture that also in higher orders they are constrained by the non-abelian exponentiation 
theorem [41, 42], as is the case for the cusp anomalous dimension. This would imply that the 
Casimir scaling relation d q JCp = d 9 JCA continues to hold at least to three-loop order. Since 
the cusp anomalous dimension obeys the same relation, Casimir scaling to three-loop order 
holds for the entire F qq and F gg functions, as shown in (17). Note that there are arguments 
indicating that for the cusp anomalous dimension Casimir scaling should hold at four loops 
and perhaps even to all orders of perturbation theory [37]. 

4 Resummation and Fourier transformation 

In the differential cross section (18) and the expression for the hard-scattering kernels C qq ^ij 
in (23) the dependence on the scales M 2 and x\ is factorized explicitly. However, for any 
given choice of the renormalization scale \x these expressions contain large logarithms, which 
need to be resummed to all orders in perturbation theory. This is accomplished by solving RG 
equations for the various component functions. We will now discuss these solutions in detail. 
We will then derive a compact, all-order formula for the hard-scattering kernels C qq ^ij in q? 
space, which is free of large logarithms. 

In practice, the easiest way to perform the resummation of large logarithms is to choose \i 
of order a typical hard-collinear scale, i.e. \x ~ qr or // ~ x^. 1 , and then evolve the hard function 
and the PDFs from appropriate initial scales to the scale /j,. The PDF evolution is standard 
and can be taken from any package that generates parton distributions. The evolution of the 
hard function reads [31] 



C v (-M 2 ,i2) = exp[2S(i2 h ,i2)-2a iq (i2 h ,i2)}i — ^- ] C v (-M 2 ,/2 h ), (52) 
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where n\ ~ — M 2 is a hard matching scale, at which the value of Cy is calculated using 
fixed-order perturbation theory. At one-loop order [30, 31] 

^(-M^^) = l + ^^(-L 2 + 3L-8 + ^ + ..., (53) 

where L = ln(— M 2 /{i\). The two-loop correction can be found in [31]. The advantages of 
using a time-like scale choice (/i 2 < 0) for time-like processes such as Drell-Yan production 
were emphasized in [51, 52, 53]. The Sudakov exponent S and the exponents a n are given by 

[54] 

S(u,n) = - — ln^rf usp (a s (/i)), a T {y,ix) = - I — rf usp (a s (/x)) , (54) 
Jv A 4 v Jv A* 

and similarly for the function a 7 ? . Note that for gluon-initiated processes the functions S and 

ar would (at least up to three-loop order) simply be rescaled by Ca/Cf- The solutions to the 

evolution equations (16) for the functions F qq and B q / N can be expressed in terms of the same 

objects. They are given in Appendix B, where we also collect the perturbative expansions of 

the anomalous dimensions and the resulting expressions for the evolution functions valid at 

next-to-leading order in RG-improved perturbation theory. 

Besides the large logarithms resummed in the evolution of the hard function |Cy| 2 , there 

are large logarithms arising from the collinear anomaly, which are contained in the first factor 

on the right-hand side of (15). In order to discuss the resummation of these logarithms in a 

systematic manner, we first rewrite 

5^) =(£) «- l ^~'- e*> 

While this formula is exact, it does not yet provide a convenient basis for a perturbative 
analysis of the resummed kernels in RG-improved perturbation theory. Counting as usual the 
large logarithm ln(M 2 //x 2 ) like l/a s , we see that the first term in the perturbative series (47) 
for the exponent F q q needs to be evaluated in exponentiated form, while from two-loop order 
on a perturbative expansion can be applied. We thus split 

F q ,(L±, *s) = ^L ± + J2 d q n (L ± ) (g) ee -&L [L ± + f qq (L ± , a.)] , (56) 

n=2 

where the function f qq starts at 0(a s ) but only receives contributions of two- loop order and 
higher. We now have 



\/r2\ -F qq (L ± ,a a ) , 2 ,,2 \ ~Vf(M 2 ,h) 

^_) _ / X T fi \ e -VF(M*,rif q9 (L ± , aa ) ? ( 57 ) 



where 



2 C F a s M 2 
r] F (M,/j,) = ln ^ (5? 

7T fi z 



counts as an 0(1) quantity as long as fi 2 ^C M 2 . 
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Our last task is to perform the Fourier transformation (23) from Xt space to transverse- 
momentum space. Since all our functions depend on xt either via a power, like in (57), or via 
the logarithm L± defined in (34), it suffices to consider the relation 

-1 [d 2 x , e— L" {^-\ ~" - i-d V - (£) " r(1 ~^ } (59) 

While gr serves as an infrared regulator, the integral converges in the ultraviolet, for xt — >■ 0, 
only as long as 77 < 1. In our formulae below we will always assume that the condition 
7]F{M 2 ,fi) < 1 is fulfilled. 4 For M = Mz or My/, for example, this implies that \x should be 
larger than about 2 GeV. Note that the higher-derivative terms in (59) are accompanied by 
powers of 1/(1 — 77), so that for 77 very close to 1 a reorganization of the perturbative series 
becomes necessary. This will be discussed in detail in a forthcoming article [44]. 

Using the above result, we obtain a closed- form expression for the resummed hard-scattering 
kernels, which reads 



C q q_> ij (z 1 ,z 2 ,qT,M 2 ,n) = \C V (-M 2 ,n)\ I q ^ i (z 1 ,-d v ,a s )Iq^ j (z 2 ,-d rv a s 

x E q q(-d v ,a s ,r]F) * " 5 



(60) 



V=Vf 



q\ \n 2 ) e 2 ™Y{r\) 

where r/p = t]f(M 2 , //), and the arguments of the Ii<-j(z, L±, a s ) functions are those shown in 
(41). It is understood that for \Cy\ 2 one uses the resummed expression in (52). All remaining 
quantities have perturbative expansions in powers of a s = a s (fj,) free of large logarithms. In 
writing the above result we have introduced the function 



E q9 (L ± , a s , 77) = exp - L ± F qq -(L ± , a s ) - 77 f q g(L±, a 
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(61) 
+ 0(a 2 ) , 



which is completely determined in terms of F qq . The two-loop coefficients d\ and Tf enter here 
already at next-to-leading order in a s . For a consistent resummation at NNLL order (or next- 
to-leading order in RG- improved perturbation theory), we need the one-loop expressions for 
the matching coefficients Cy and ij«-j, the two- loop expression for the exponent F qq , the two- 
loop expression for the anomalous dimension / ~f g , and the three- loop cusp anomalous dimension 
and /3 function. All of these ingredients are known. 

Note also that, owing to the simple q\ dependence of the resummed result (60), it is trivial 
to perform the integral over transverse momentum required to calculate the cross section (27) 
defined with a cut on transverse momentum. 

5 Asymptotic divergence and reorganized expansion 

Despite the fact that it correctly resums all large logarithmic terms in the perturbative series, 
the elegant formula (60) just derived is of limited practical use. The reason is a strong factorial 



4 The singularity at 77 = 1 was noted a long time ago and referred to as a "geometrical singularity" [55]. We 
will show in Section 5 that it disappears when a class of factorially divergent higher-order terms is resummed. 
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divergence of the perturbative expansion coefficients resulting from terms in the functions I^j 
and E q g of order (a s L^) . To understand the origin of this effect, we recall from (61) that 
before expansion in powers of a s the hard-scattering kernels contain quadratic terms in L± 
in the exponent. The same is true for the C?[(a: s LjJ ] terms in the kernels U^j, which 
exponentiate as a consequence of the cusp logarithm in the evolution equation (22). Let us 
then consider, instead of (59), the Fourier integral 



2- ke-*^e*44 = £11 r M jJe^b^)e^-\^ = £ll K L a to 
Air J /i 2 J_ OQ \ ix f fi 2 \ fi 2 

(62) 

where bo = 2e 7B , and in the case at hand 



2n 



[T^ + VF (M 2 ,fi)f3 ]. (63) 



Some useful properties of the function K(r], a, r) are summarized in Appendix C. The above 
definition is such that for a = we recover, up to a trivial factor, the result (59) with n = 0: 

Kin r) = r^ 1 "~ r} > — (64} 

^VW) e 2(ij-ihj5 r(rj) ' 

Keeping the quadratic term in the exponent vastly improves the convergence behavior of 
the Fourier integral. For a = (i.e., without the quadratic term) the integral on the left-hand 
side of (62) converges in the ultraviolet (for xt — > 0) only if rj < 1, and for rj < | its value 
must be defined by analytic continuation. For a > 0, on the other hand, the integral converges 
for all values of r). It is then perhaps not surprising that any attempt to expand the Gaussian 
weight factor in a perturbative series leads to a badly behaved expansion. Indeed, writing the 
formal series 

*(w)L, = t h (-;)"fr*(i.M = E k (-i)> r "' *ff A ' (65) 

n=0 n=0 V " 

it is not difficult to see that the series is factorially divergent. To illustrate this point, we 
consider the special case where r = 1 (corresponding to the default scale choice \i = qx) and 
rj is close to the critical value 1. One then has 

r(1 -' ) - ' 2 S-^-^(l -,)* + •.., (66) 



e 2(„-i) 7B r(?7) I-77 3 v " 5 

and taking 2n derivatives of the leading term generates (2n)!/(l — rj) 2n+l . A more careful 
analysis reveals that 
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+ J2k n a n + 0(l- V ), (67) 

?1=0 



where the coefficients k n do not exhibit the strong factorial growth of the terms in the first 
sum. While this series is badly divergent, the fact that it has alternating sign implies that it 
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can be Borel-summed. We obtain 
K(rj,a,l) 
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n=0 



where Erf(x) is the error function. Note that the singularity at r\ — 1 has disappeared after 
Borel summation. Expressions for the first few k n coefficients can be readily derived in terms of 
7s and ( n values. Numerically, we find k « 0.3152, k x « 0.2431, A; 2 « -0.0747, A; 3 « -0.0509, 
/C4 ~ 0.0446, /C5 ~ —0.0043, .... We have checked numerically that the Borel-summed result 
(68) reproduces the original integral as defined in (62) precisely. 

The presence of factorially divergent terms in the perturbative expansion of the CSS for- 
mula for the resummed Drell-Yan cross section at small transverse momentum was first pointed 
out in [55]. These authors considered this a problematic feature of qr resummation, which 
arises because the integral over xt in (23) includes the region xt < M~ 1 , where the resumma- 
tion formula is not justified. For the same reason, the authors of [21] proposed a modification 

2 2 2 

of the formula in which one replaces In ^^ > l n [(|tf + apO/^ 2 ], thereby effectively introduc- 
ing a UV cutoff on the xt integral. Note, however, that the quadratic term in the exponent 
of the Fourier integral (62) provides a Gaussian cutoff ensuring that \L±\ < l/y/o,, which 
automatically excludes values xt < M -1 . 5 Our interpretation of the origin of the divergent 
terms is a different one. The asymptotic behavior of the perturbative series in gauge theories 
is a hint that perturbation theory by itself can only provide an incomplete description of na- 
ture. For instance, the appearance of infrared renormalon singularities - factorially divergent 
terms related to running-coupling effects, whose Borel sum is ambiguous by power-suppressed 
amounts - hints at the existence of non-perturbative power corrections, which must be added 
to the perturbative series in order to obtain a precise description of strong-interaction physics 
[56]. In the present case, the appearance of a factorially divergent series indicates that some 
essential features of the transverse-momentum distribution in the region of small to moderate 
q? cannot be reproduced at any fixed order of perturbation theory, because they carry some 
genuinely non-perturbative information (even though this does not necessarily indicate a sen- 
sitivity to long-distance physics). It is well-known that for sufficiently large M the intercept of 
the distribution da jdq\ is calculable and depends in an intrinsically non-perturbative way on 
the coupling constant, da/dq^\ qT=0 oc e const//as /\/^I [2, 4]. Our Borel-summed expression (68) 
shares many of the features of the formula for the intercept. Only after Borel summation of 
the divergent series, which is possible due to its alternating sign behavior, we obtain a reliable 
description of the qx spectrum, which encodes its non-perturbative behavior in an appropriate 
way [44]. 

Let us now discuss the implications of these findings for the resummed Drell-Yan cross 
section. The main lesson is that in order to avoid an asymptotic growth of the perturbative 
expansion coefficients one should keep the quadratic terms in L± in the exponent of the Fourier 



5 For 77 < 1 the integral is UV convergent even without the quadratic term, and contributions from the 
region ly < q^, are power suppressed. 
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integral (62). This gives rise to a new formula for the hard-scattering kernels in (23). Instead 
of the expression (60), we obtain 



C q q^ ij (z 1 ,z 2 ,qT,M 2 ,ij,) = |CV(-M 2 ,/i)| I q ^ i (z 1 ,-d r) ,a s )Iq^ j (z 2 ,-d ri 
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r 1 =r lF (M 2 , l j,) 



(69) 



where a is given by (63), and the functions I^j and E qq are obtained simply by dropping the 
0(a s L 2 ± ) terms in the expressions for ij_^- and E qq in (41) and (61), and likewise in higher 
orders. The integral K(r), a, r) and its derivative with respect to rj can easily be evaluated 
numerically. We have not succeeded to derive a suitable analytic expression for this integral 
in the general case where r ^1 and rj is not close to 1. 



6 Comparison with the literature 

The standard formalism for transverse-momentum resummation has been developed in a sem- 
inal paper by Collins, Soper, and Sterman (CSS) [4]. According to this work, the resummed 
differential cross section at leading power can be written in the form 
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where \x b = b^/xx is assumed to be in the perturbative domain. It is a straightforward exercise 
to work out the relations between the various objects in this formula and ours. We find 
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The one-loop coefficients are d\ — and 



e? = Cf ( ^ - 16 ) . (73) 



The two-loop coefficient g?2 has been given in (51), while e\ can be extracted from the results 
compiled in [31], however it contributes to B(a s ) at 0(afy only. We have checked that the 
relations in (71) are compatible with our perturbative results. 

Note that according to (71) the coefficient A in the CSS formula differs from the cusp 
anomalous dimension starting at three-loop order, and the coefficient B differs from the quark 
anomalous dimension 2 r y q starting at two-loop order. 6 The first non-zero deviations are (here 
A^ and B^ denote the n-th order coefficients in the expansion in powers of a s /(47r)) 

A< 8 > = if + 2fa<% , B^ = 2 7 ? + d\ + /3 e? . (74) 

The two-loop expression for B(a s ) was obtained a long time ago in [6], while for gluon-initiated 
processes such as Higgs-boson production the corresponding coefficient was calculated in [7]. 
Using these results, we have derived the anticipated relation (51). Inserting the coefficients 
d^' 9 into (74), we obtain the coefficient A^ 3 \ which up to now was the last missing ingredient 
for a full NNLL resummation of the q? spectrum. In the literature it is often assumed that 
A^ = rf (see e. g. [21, 22]), which is true for soft gluon resummation, but our results 
show that for transverse-momentum resummation an extra contribution arises because of the 
collinear anomaly. Numerically, for the quark case with n/ = 5, we find rf = 239.2 while 
A^ = —413.7, so the extra term is much larger than the contribution from the cusp anomalous 
dimension and has opposite sign. It will be interesting to see how this changes the numerical 
predictions for the spectrum. Note also that, due to Casimir scaling, in the gluon case a similar 
situation but with larger coefficients occurs, and we find T^ = 538.2 while A^ = —930.8. 

CSS also derive a formula for the resummed cross section at small qx in terms of transverse- 
position dependent PDFs Pj/jv^^t), which is equivalent to our result (18) once we identify 

V i/N (i } x T ) = \C v (-/j, 2 b ,/j, b )\B i/N (£,xl,/j, b ). (75) 

The definition of the function Vi/N^, %t) proceeds in two steps. One first solves a differ- 
ential equation governing the ( dependence of the gauge-dependent distribution function 
Vi/N(£,XT, A*; C) mentioned earlier in our discussion in Section 2.2. This introduces a gauge- 
independent boundary function V(£,XT,fj) [34]. The function Pi/Ar(£,a>r) is obtained by 
multiplying this function with an a^-dependent factor that cancels its scale dependence [57]. 
Like in our case, no operator definition for the final transverse PDF Vi/n is provided. 

While we obtained the closed- form expressions (52), (60), and (69) for the cross section 
directly in momentum space, the CSS formula (70) involves a Fourier integral over xt- The in- 
herent scale choice \x-\X\, — bo/xr eliminates all L± logarithms and thus automatically resums 
the factorially divergent terms discussed in Section 5. However, since the integrand involves 
a s (/ib), the integration hits the Landau pole in the running coupling, so that a prescription is 



6 The first relation in (71) can be found, in almost precisely this form, in equation (3.13) of [4], from which 
it follows that F q q(xj-,fi) = —K(xT^,a s ) in the notation of that paper. 
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needed to regularize the integral. In practical applications, the integration is cut off at large 
xt values. To account for the missing long-distance contributions a non-perturbative model 
function is used. For q? in the perturbative domain these contributions are formally power 
suppressed, but it is irritating that an explicit prescription for how to deal with them is needed 
even for q? values deep in the perturbative regime. Special care has to be taken in order not 
to induce unphysical power corrections in the cut-off procedure [58]. The explicit cut-off also 
makes it difficult to perform the matching to fixed-order computations, since resummation 
effects persist even when the CSS formula is evaluated at large qx ~ M, as discussed in [59]. 
All of these complications are absent in our resummed results (52), (60), and (69). Our ex- 
pressions are given directly in momentum space and do not involve a Landau-pole singularity. 
In the spirit of effective field theory, we never perform scale setting inside integrals over the 
running coupling a s (fi). Instead, the scales are chosen such that the integrated result is free of 
large logarithms. The matching onto fixed-order computations is completely trivial, since our 
analytic result (60) can easily be reexpanded in powers of a fixed coupling a s (/x). Moreover, 
in contrast to the CSS formula, the resummation switches itself off adiabatically when /i ~ q? 
approaches M, since all the logarithms become small and the RG evolution stops when these 
scales become equal. 

The advantages of performing the resummation in qx rather than in xt space were empha- 
sized a long time ago in [59], where an improved leading- log formula extending the work of 
[1] was derived, which resums all logarithmically-enhanced terms to two-loop order. We have 
checked that this formula is consistent with our exact result (60) provided the two-loop mod- 
ification of the B^ shown in equation (34) of [59] is included. Note that the approximation 
adopted in that paper misses terms of order (a s L) n L™~ 4 for n > 3 (denoting L = \n(M 2 / q 2 -,)) , 
which for a s L = 0(1) become progressively enhanced for increasing n. In our approach, all 
such enhanced terms are resummed. 

Transverse-momentum resummation has been studied earlier using SCET in [23, 24] and 
more recently in [25]. These papers solve the RG evolution equation of the hard function \Cy\ 2 
in a way that is equivalent to (52). The exponentiation of the additional large logarithms of 
q^/M 2 , which arise from the transverse PDFs due to the collinear anomaly, is however not 
addressed in any of these works. Because of this, the relations derived in [23, 24] between 
the effective-theory results and the traditional CSS approach miss the extra contributions 
proportional to g\ in (71) and hence are only valid to next-to-leading logarithmic order. In 
[25] the resummation of large logarithms was studied in momentum space, similar to our 
approach. However, since their result is given in terms of a multi-dimensional convolution 
integral it is more complicated to implement in practice. More importantly, the fact that 
the large logarithms arising from the collinear anomaly (which would have to appear in the 
matching conditions for the jet and soft functions in that paper or in the convolutions of 
these objects) were neither addressed nor resummed implies that the resummation formulae 
derived in [25] will have to be modified already at NLL order (with the usual counting of 
logarithms in the exponent). The reason is that they do not accomplish the resummation of 
the 7^- dependent terms in (60). 

At first sight, the papers [23, 24, 25] avoid dealing with the collinear anomaly by keeping 
power-suppressed terms in the computation of the soft function and of the matching of the 
transverse onto the usual PDFs. In particular, the multipole expansion of the soft function 
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is not performed, and for this reason the soft function (8) survives in the final factorization 
formula. While keeping higher-order terms regulates the light-cone singularities, care must 
to be taken to avoid double counting (this requires so-called zero-bin subtractions). Also, if 
the multipole expansion is not performed consistently, the soft and beam functions depend on 
several scales, such that their computation is not only much more cumbersome than in our 
case, but also gives rise to extra logarithms, which would need to be resummed. No attempt 
is made in [23, 24, 25] to derive the corresponding evolution equations. 

7 Conclusions 

We have analyzed the production of high-mass Drell-Yan pairs at low transverse momentum q? 
using effective field-theory methods. Since the classical paper of Collins, Soper, and Sterman 
from 1985, it is known how to resum the large perturbative logarithms of qr over the invariant 
mass of the lepton pair which arise in this kinematical region. Despite this, certain aspects 
of the factorization properties of the cross section at small qx have continued to be puzzling. 
The difficulties arise because the cross section factors into a product of transverse-position 
dependent PDFs, whose naive definition leads to inconsistencies. Many different "improved" 
definitions for these objects were proposed over the past thirty years, but none of them seems 
entirely satisfactory. 

In this paper we have shown that while individual £<r-dependent PDFs are not well defined 
without additional regulators, their product is. However, this product depends on the large 
momentum transfer q 2 ^> l/x T of the underlying hard-scattering process. This anomaly arises 
because the regulators needed to define the individual PDFs in SCET necessarily break the 
naive factorization of the effective Lagrangian into two independent collinear sectors. Similar 
to other anomalies, the breaking of factorization has very specific properties, which imply 
that the dependence on the momentum transfer exponentiates. Only after the q 2 dependence 
associated with the anomaly is accounted for, the remainder factors into a product of two 
functions, each dependent on the transverse separation and the fraction of the large momentum 
carried by the leading parton. These functions provide an operational definition (though not an 
operator definition) of transverse-position dependent PDFs, whose evolution equation is solved 
exactly and is known explicitly at three-loop order. However, it is an open question whether 
the anomalous q 2 dependence is truly universal and process independent. We note that at 
large xt ~ 1/Aqcd the exponent of the q 2 dependence from the collinear anomaly becomes a 
non-perturbative quantity, so that the dependence on the large momentum transfer q 2 ^> Aq CD 
is no longer perturbatively calculable. For small transverse separation iy <C 1/Aqcdj on the 
other hand, the transverse PDFs can be matched onto standard PDFs, and their dependences 
on xt and q 2 are calculable in perturbation theory. We have performed the necessary matching 
calculation at one-loop order and derived the q 2 dependence due to the anomaly at the two- 
loop level. 

As a result of our factorization analysis, we are able to obtain the first, closed analytic 
expression for the cross section at small q? <C M directly in momentum space, which is free 
from large perturbative logarithms. In addition to logarithms arising from the evolution of the 
hard function to lower scales, which we resum by RG evolution and which also arise in soft- 
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gluon resummation, our result contains a second source of large logarithms due to the collinear 
anomaly. A crucial result of our analysis is that these additional logarithms exponentiate. We 
have also shown that, besides parametrically large logarithms, the perturbative series for 
the transverse-momentum distribution contains terms exhibiting a strong factorial growth at 
higher orders in a s , which must be summed in order to obtain reliable results. We have 
explained the origin of these terms and shown how they can be accounted for to all orders in 
perturbation theory. 

Our momentum-space expression for the resummed cross section offers a number of ad- 
vantages over the traditional formalism. Since we have an analytic result for the resummed 
cross section, we can easily reexpand it to match to fixed-order results, while the matching 
is nontrivial in the traditional approach. Our resummation automatically switches off as q\ 
approaches q 2 . In addition, our result is free from unphysical Landau-pole singularities, so 
that no explicit prescription for the non-perturbative regime is needed as long as qr ~^> Aqcd- 
We have derived all the necessary input to perform resummation at next-to-next-to-leading 
logarithmic accuracy and as a by-product have obtained the three-loop coefficient A^ 3 \ which 
is needed in the traditional approach to achieve this accuracy. A phenomenological analysis 
of electroweak boson production using our formalism will be presented elsewhere. 
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A Generalization to W and Z production 

The generalization of our results from the photon-induced Drell-Yan process discussed in the 
paper to W and Z production is straightforward. To obtain the double differential cross 
section d 2 a jdq\dy from (18) and (24), one needs to change the prefactor according to 
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and insert the proper charge factors. For the Z boson one must replace 
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where 6w is the weak mixing angle. Since the W bosons have flavor- changing couplings, to 
sum over flavors must be replaced by a double sum over individual quark and anti-quark 
flavors, q and q' . Only left-handed currents appear in this case. The relevant coupling for a 
W~ boson produced in the annihilation of an anti-up and a down quark is 
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where V q > q are elements of the quark mixing matrix. 

B RG evolution 

Here we give the perturbative expansions of the functions S and ap defined in (54), working 
consistently at next-to-leading order in RG-improved perturbation theory. At this order we 
need to keep terms through 0(a s ) in the final expressions. The result for ar is given by [54] 
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A similar expressions, with the coefficients Tf replaced by 7^, holds for the function a 7 <? 
expression for the Sudakov exponent S is more complicated. It reads [54] 
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where r = a s ([A)/a s {v). Whereas the two- loop anomalous dimensions and /3- function are 
required in (Bl), the expression for S also involves the three-loop coefficients rf and /?2- 

The relevant expansion coefficients of the cusp anomalous dimension r^ sp to three-loop 
order are Tq = ACp and [60] 

rf /67 tt 2 \ 20 
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Due to Casimir scaling, these ratios are the same in any representation of the gauge group. 
The coefficients of the quark anomalous dimension 7 9 to two-loop order read [31, 37] 

H = -3C F , (B4) 
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The three-loop coefficient 7I can be found in the same papers. Finally, the expansion coeffi- 
cients for the QCD /3-function to three-loop order are 

Po = y c a - g 7>n/ , 
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We finally present the exact solutions to the RG equations for the exponent F qq - and the 
function B q / N given in (16). The solution to the first equation reads 

F Q q(xT, /i) = -2a r (/i , y) + F qq (xT, /io) = -2a r (// b , //) + J^ < ( ^ ) , ( B 6) 

n=l ^ ' 

where /io is an arbitrary initial scale of order qx or x^ 1 , and in the last step we have made the 
special choice /x = Hb — 2e~ TB /a>r. The solution to the second equation is 

/ „2 2 \ -2ar(M0,M) 

B i/N (z,x%,(i) = exp [-25 , ( < u ,/i) + 2a 7g (/i ,/i)] ( ^^ j B i/N (z,x^, n ) . (B7) 

As before, a particularly convenient choice is to set /io = A^. This leads to 

B q/N (z, 4, /i) = e ^(4.M) 5^^ 4) , (B8) 
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where h qq (x^,ii) = — 2£'(/x&, //) + 2a 7 ? (//&, /i) , while B q / N (z,x^) = B q / N (z,x^, Hb) is RG invari- 
ant. This is the function that, up to a factor, is equal to the distribution function V q /N(z, %t) 
introduced in [57], see (75). For /i ~ fi^ one can construct a fixed-order expression for the 
exponent h qq as a polynomial in L±. At two-loop order we obtain 
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C Properties of the function K(rj,a,r) 

Explicit expressions for the Fourier integral can be obtained in the limits of large a or small 
r. We obtain 
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a> 1, 



and 



K(r},a,r) 




An 
a 



- e 1E re a + C(r ) 



e a 



r < 1. 



(CI) 



(C2) 



Noting that a rescaling of the ratio qr/fJ' in the argument of the Bessel function in (62) 
can be compensated by a shift of the integration variable £, we find that 



K(n, a, r) = r ^' 1 ~-* lnr K (n - | In r, a V 

where K(s, a) = i^(s, a, 1). This function obeys the partial differential equation 

d 2 d \ - - T(l - s) 

— + 4—)K(s,a) = 0, K(s,0) = ^r r \ — -y- . 
ds 2 da) y ' ' v ' ; e 2 ( s ~ 1 )^r(s) 

The Borel image of the divergent part of the series (the first sum) in (67) is 



(C3) 



(C4) 



BmJu) 



l — rj 



S(u)- 



1 — 7] 



2( u + (i- v yf 2 



-2 7 B 



8(u)- 



1 



2 (1 + uf 2 



(C5) 



and performing the Borel integral f^ du B^ v (u) e~ u l a yields the expression shown in (68). 
Note also that relation (CI) implies that for large a the non-divergent terms in the series must 
obey the sum rule 



7 , k r . 

n=0 



2J (6 ) - 1 + e' 2TE +0(a 



,-3/2> 



1.2988 



(C6) 
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